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We calculate the complete double logarithmic contribution to cross sections for semi-inclusive 
hadron production in the modified minimal-subtraction (MS) scheme by applying dimensional reg- 
ularization to the double logarithm approximation. The full double logarithmic contribution to the 
coefficient functions for inclusive hadron production in e + e~ annihilation is obtained in this scheme 
for the first time. Our result agrees with all fixed order results in the literature, which extend to 
next-next-to-leading order. 
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CN ; I. INTRODUCTION 

. The inclusive production of particles in the framework of the factorization theorem and perturbation theory at 
high energy has been understood since a long time. However, perturbation theory fails when the fraction x of 
O | available energy carried away by the observed particle is too low. Specifically, large unresummed logarithms spoil the 

■ convergence of the series. The wealth of data at lower x that have to be excluded from global fits pending the explicit 
' resummation of such logarithms would significantly improve the constraints on fragmentation functions (FFs) at large 

x and, for the first time, at small x, as well as on a s (Mz). 
^vq ' In fact, the largest logarithms, the double logarithms (DLs), in the splitting functions that determine the evolution 
, of the FFs have been determined to all orders a long time ago [l[ , and have even been used to perform leading order 
(LO) global fits to data measured at the smallest x values. Specifically, in these fits the evolution was calculated 
in the fixed order (FO) approach, to allow for a description of the large x data, while including the complete DL 
contribution to all orders, and this consistent approach resulted in a simultaneously good description of the remaining 
smaller x data. 

Strictly speaking, the theoretical approach used at LO in the global fits of Refs. 0, Q to measurements of inclusive 
particle production in e + e~ collisions is incomplete, because the DLs in the coefficient functions are not resummed. 
These DLs are expected to be not as important as those appearing in the evolution because they only appear at 
and beyond NLO and, furthermore, the approach turned out to be adequate for the numerical analysis. However, 
the inclusion of the DLs in the coefficient functions could make a significant improvement to the accuracy of cross 
section calculations making the analysis of Refs. 0, H[ feasible also at next-to- leading order (NLO). The complete 
DL contribution to partonic cross sections has been calculated [3] for the case in which the collinear singularities 

■ are regularized by giving a small mass m g to the gluon, the so-called massive gluon (MG) regularization scheme. 
This result turns out to be inconsistent with the full next-next-to-lcading order (NNLO) result in the modified 
minimal-subtraction (MS) scheme. In the literature, the 0(a 2 s ) contributions to the fragmentation functions in e + e~ 
annihilation have been originally computed in Refs. 0-0] • Their results were confirmed 10 years later using rather 
different technologies by the new independent computation in Mellin (N) space Q. The Mellin space result of Ref. 
§ agrees with the Mellin transform of the x-space result of Refs. 0-0|_performed in @. As a result of these last 
computations, several typographical errors present in the original papers |5|-|7[ have been corrected. The inconsistency 
between the NNLO DLs calculated in Ref. Q and those calculated in Refs. 0-0, E3] is not surprising — the two 
computations were carried out in two different regularization and factorization schemes. 

Therefore, it is necessary to calculate the DL contribution to cross sections also in the MS scheme, which is the 
goal of this paper. In section UH we review the derivation of the double logarithm approximation (DLA) result 
calculated in the MG scheme, and its factorization, for the probability of the inclusive production of an observed 
gluon. Indeed, it is the gluon channel which contains the DLs. In section UlTl we calculate the modification to this 
result when dimensional regularization (DR) is used instead, and factorize it in the MS factorization scheme to obtain 
the complete DL contribution to the splitting functions and coefficient functions in this scheme. Finally, in section 
W] we present our conclusions. 
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II. DOUBLE LOGARITHMS IN THE MASSIVE GLUON SCHEME 



We begin by repeating the derivation of the result for the DL contributions to splitting functions and coefficient 
functions of Ref. [J] in a way which makes suitable its comparison with our derivation using dimcnsionally rcgular- 
ization, to be presented in section IIII1 In this section we work in the MG scheme in order to regularize the mass 
singularities, which is why we assign a small mass m g to each gluon. 

As is well known, these DL contributions appear in the gluon-gluon and gluon-quark timelike splitting functions 
[l[ and in the timelike gluon coefficient function To extract them we consider a general process with a colour 
singlet final state involving the production of an "observed" gluon of momentum q from a hard parton of momentum 
p around which a jet is formed. In the DLA, the DL contribution arises from unobserved soft gluons in the final state. 
Consequently, there must be an additional hard parton to account for the recoil from the parton of momentum p as a 
result of momentum conservation. The momentum of this additional parton will be written p. The cross section for 
this process will be written da(p,p,q). A typical example is the process e + + e~ — > V* — > Q(p) + Q(p) + g(q) + X, 
where V* = 7, Z is a virtual vector boson, where the jet is formed around the quark Q with momentum p and around 
the antiquark Q with momentum p, and where X is any hadronic final state that is allowed by quantum number 
conservation. Thus, to obtain the DL contribution to the cross section, we consider the configuration in which the 
unobserved part consists of only N soft gluons of momenta qi, q 2 , • ■ • , In whose phase space is fully integrated out. 
Therefore, defining da]\[(j),p, k\, k 2 , ■ ■ • , fc/v) to be the cross section in which N gluons of momenta k a , a = 1, 2, . . . ,N, 
are produced together with the partons of momentum p and p, we can write 

00 

d<r(p,P,q) = ^2 da N+1 (p,p,q,q 1 ,q 2 , . . . ,q N ), (1) 
n=o 

where it is understood that the q a are fully integrated over, but not q. It is a well known result fill ] that the DL 
contributions come from the kinematic configuration in which the momenta of the soft gluons are strongly ordered, 
i.e. 

\q\ « |«!| « |<f 2 | « • • • « \q N \ « Q/2, (2) 

where Q 2 is the perturbative scale, i.e. the scale that the factorization and renormalization scales should have the same 
order of magnitude as. For a general process, the choice Q 2 = (p + p) 2 is suitable and is usually made, in particular 
this choice is made in e + e~ annihilation. Furthermore, as initially explicitly proved up to three loops [l| and then 
formally at all orders in Ref. [ll|, to extract the most singular terms the phase space should also be restricted to the 
region where the angles 0/ of the emitted gluons with respect to the hard parton of momentum p are strongly ordered, 
i.e. 



9 < 0i < 2 < • • • < 9 N < 1, (3) 

where refers to the gluon of momentum q. Now, in Appendix \K\ we derive the factorization of soft gluon emission 
without carrying out the phase space integrations for any of the soft gluons. Thus, using Eq. (|Al3| . which is the 
general result for the factorization of soft gluon emission, but choosing N + 1 final state gluons and setting k\ = q 
and k a+ i = q a for a. = 1, 2, . . . , N, da^ + i at DL accuracy factorizes in the following way: 

da N+ i(p,p,q,q 1 ,q 2 ,...,q N ) = da B (p,p)dWg(q)dwg(qi)dw g (q 2 ) . . . dw g (q N -i)dw I (q N ), (4) 

where dasiPiP) is the Born cross section of the underlying process — e.g. in e + e~ annihilation above it is the LO 
cross section for the process e + e~ — > Q(p) + Q(p), and I labels the species of the initial state parton with momenta 
p (see Appendix According to Eq. (|A10j) . dwi(q a ) is, in d = 4 dimensions, given by: 

dwiiqa) = 9 Kj iv- qa ){P- qa )wm^ (5) 

where g is the strong coupling constant and Kj is defined in Appendix [X] We will work in a frame where p = — p. 
Aligning the z-axis with p, this means that 

p=|(l,0 ± ,l), p=|(l,0 ± ,-l). (6) 

We introduce the usual Sudakov paramctrization for the gluon momenta q a , i.e. 

q a = x a p + z a x a p + q a ±, q a ± = (0, q Q j_, 0), q 2 a± = z a x 2 a Q 2 - m 2 , (7) 
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where we have used the on-shcll condition for the gluon in the MG regularization scheme (i.e. q\ — m 2 ). From Eq. 
([7]), we find that the quantities x a and z a are given by 

_ v q* _ i - \q a \/q a cos6» Q 



p-q a 1 + \q a \/q° cost 
p-p Q(l + Z a )' 



(9) 



9° = Vl<7«l 2 +m2. (10 ) 

We can now use Eqs. ([5} and © to change variables in Eq. (|5|) from g Q , p and p to x Q and z a . Calculating the 
Jacobian determinant of this change of variables and using d 3 q a = 2ir\q a \ 2 d\q a \d cos 9 , where the factor 2ir arises from 
the azimuthal integration using the symmetry of dwi around the z-axis, the result is that 

1 / \ T jr ds *X> d' 2s (~y _ _ . Q , , v 

dw I {x a , z a ) = 2a s Kj h 0{m g ), (11) 

where a s = a s /(27r) = g 2 /8ir 2 is the bare coupling, which at this level of accuracy can be replaced by the renormalized 
coupling. In Eq. (fTTj) . we have neglected terms proportional to m 2 which do not affect the DLs we are interested 
in. However, in the MG regularization scheme, the small gluon mass m g results in a non-zero lower cutoff for the 
variables z a appearing in the denominator of Eq. (jll[) . To find the value of this small cutoff Z a , we note from Eq. 
([8]) that the minimum of z a corresponds to the collinear limit, i.e. a = 0. In this limit, using Eq. ([9]) together with 
the on-shell condition in Eq. (|10p to express in terms of x,m g and z, Eq. ([5]) leads to the following implicit 

equation for the cutoff Z a : 



This has the exact solution 



1 J-\ 4m 3 

(12) 

1+ 4 /1 4m2 9 



Z « = -T7»- ( 13 ) 



The physical meaning of the variables x a and z a defined in Eqs. and (J9j) becomes clearer if we keep the first order 
term in the small gluon mass m g and in the small angles a as required by Eq. ([3]): Defining tu a = q&, 



Po 





(14) 



It is clear that if we use these expressions to define x a and z a appearing in Eq. (|11[) , together with the collinear cutoff 
in Eq. ((13)) to compute the DLs according to Eqs. (U} and ((H), only 0(m 2 ) terms are modified, but not the terms 
which are singular or just non-vanishing as m g — > 0. We may therefore set xq = x, where x is the usual energy or 
momentum fraction, zq = z and Zq = Z . Since x is measured in experiment but the remaining degrees of freedom in 
q are not (i.e. z is integrated over), we replace dw g (x, z) with dw g (x' , z)S(x — x')dx on the right hand side of Eq. (fTj) 
and then divide both sides by dx, giving 

d 

— (p,p,x, Z) = da B {p,p) ^2 S(x' - x)dw I (x',Zo)dWg(x 1 ,zi)dWg(x2,z 2 ) ■ ■ ■ dw g (x N , z N ). (15) 
x n=o 

In Eq. (|15[) . we have noted that there is an explicit dependence on m g , which determines the lower bound in the 
integrations over z a according to Eq. (|13[) . 

Finally, putting Eq. (fTTj) into Eq. (|T5|) . then using Eqs. ^ and © to explicitly determine the integration limits in 
terms of various variables defined above including the collinear cutoff defined in Eq. (|13p , we find that 

r 1 dx 1 f 1 dz' 

xg(x, Z) = 2a s Kj I / —x\G{x 1 ,z), (16) 

Jx x i Jz z 
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where 

5( X ' Z ) = 3 — ]~^:^r(P,P,x,Z) (17) 
a<7B(P,P) dx 

is the process-dependent (i.e. /-dependent) probability density in x for the inclusive emission from an initial state 
parton of species / of an observed gluon carrying away a momentum fraction x, and 

xG(x,Z)=S(l-x)+f(2a s C A r /* — f — - t — t ~ C ~ • • • t ~ ( 18 ) 

£Ti Jx x 2 J X2 X 3 J XN _ 1 X N J z Zt J Z1 Z 2 J ZN _, Z N 

is the associated process-independent probability for the inclusive emission from an initial state single gluon of an 
observed gluon carrying away a momentum fraction x and emitted at an angle v obeying Z = (1 — cos^)/2. Notice 
that Eq. (|18p explicitly contains the conditions z a > z a _\ for a = 2, 3, . . . , N dictated by the angular ordering in 
Eq. ([3]) plus the condition z\ > Z. However, it is also compatible with the less restrictive conditions z a > Z a for 
a = 1, 2, . . . , N, as may be understood by also exploiting the conditions Z a -\ > Z a for a = 2, 3, . . . , N dictated by the 
momentum ordering in Eq. In fact, we have z\ > Z > Z\ and z a > z a _i > Z a -\ > Z a for a = 2, 3, . . . , N, so that 
it follows by iteration that z a > Z a is satisfied for all a = 1, 2, . . . , N. We note that, in Eq. (fl"6|) . z' = m 2 g / (x' 2 Q 12 ) 
plays the role of the collinear cutoff in the calculation of G in Eq. (fT5|) . where Q' 2 = zx 2 Q 2 is the energy scale relevant 
to G in the same way as Q is the energy scale relevant to g because, using Eq. ([7]), we see that Q' 2 is the transverse 
energy of the first gluon emitted from the initial state parton of momentum p. 

In order to factorize <?, we will first consider the factorization of G, which turns out to be simpler and more 
fundamental. The quantity g, which can be more directly related to physical observables (see later), can then be 
obtained from G via 

g(x,Z) = ^[G(x,Z)-S(l-x)}. (19) 

We observe that this equation, together with the definition of g(x,Z) in Eq. (|T5|. represents the integral master 
equation for G(x, Z) originally derived in the framework of a pure gluonic theory in Ref. whose solution is actually 
given by Eq. (fT8|) . 

Performing the integrals in Eq. (|T8| gives 

1 °° 1 111 

G (x, Z) = 5(1 -*) + -£ (^Aa s ) N —— ln^ 1 -- ln N -. (20) 

N=l ^ >' 

We define 



G\x,-^j=G{x,Z(x)). (21) 
In Mellin space, defined for any function f(x) via 

f(oj)= [ dxx"f(x), (22) 

Eq. (HOD reads 



o 



° (»■ # J - £< 2C ^->" g„ (N - l)'.mi.{N - rny.u N + m " ' < 23) 



where we have used Z = m 2 /x 2 Q 2 and the fact that 



1 f 1 In 1 '" 1 ! 1 

dra" = — . (24) 



(r — 1)! J a; w r 

Now, according to the QCD factorization theorem [l2|, [l^, we know that all the collinear singularities in Eq 
(which here appear as logarithms of the regulator m g ) can be factorized into a transition function T according to 

Gk^l =c( U ,a.,^)r(w,a,&), (25) 
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where 



u,a s , ^| ) 



rrig 



(26) 



with fip being an arbitrary factorization scale. Choosing \j? F = Q 2 , writing C{uj,a s , 1) = C(u,a s ) and expanding Eq. 

»2) 



5]) in powers of ln(Q 2 /m 2 ) gives 



G I u, ^§ I = C7(w, o a ) - 7 (w, oOCfw, a.) In + 0(]n'(m;/Q )) 



(27) 



Therefore, comparing Eq. (|27[) with the part of the summation in Eq. (|23[) for which m — N (i.e. the coefficient of 
In gf ) and then for which m = N — 1 (i.e. the coefficient of In 1 gf), we easily obtain that 



N=0 



(N — l)\N\co 2N 



and 



< w ^ VVor , j v (-2) JV - 1 (27V-2)! 

JV=0 



[(TV - l)!] 2 u; 2Ar 

These sums may be respectively identified to obtain, for the coefficient function, 



iw + Vw 2 + lec^os 

C(w,a s ) = -- 



2 Vw 2 + 16CUa s ' 
as well as the result determining the anomalous dimension, 

2Cao>s 



j(u,a s )C(u),a s ) 



\/u 2 + 16Cao- s ' 



(28) 



(29) 



(30) 



(31) 



7(w,a s ) = i ^-cj + \/uj 2 + 16CUa s ) . 



The quantity T obeys the timelike DGLAP equation, 



2 -^ r ' 

dm n F 



(32) 



(33) 



so that the dependence of the coefficient function on the factorization scale is easily obtained by imposing the inde- 
pendence on Hp of the gluon density function G in Eq. (f2"5j) at the same level of accuracy. 
Next we turn our attention to the factorization of g. Defining 



(34) 



Eq. (jT9|) in Mellin space becomes 



g ro 2 " 



Kj_ 

C A 



Q 2 



- 1 



(35) 



The general factorization formula reads 



(f,g) = (C q ,C g )T, 



(36) 
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where / describes the equivalent partonic process to g but for the inclusive production of a quark instead of a gluon, T 
is a 2x2 matrix containing all collinear singularities in (Q,g), and C q and C g are respectively the collinear-singularity- 
frcc quark and gluon coefficient functions. In the DLA, / = C q — 1 so that, from Eqs. ([25]) and (|36|) . we find that 



r ( Q2 \ Kl 
C„ | u,a s , — I = 



C A 



C { u,a s , — ] - 1 
Mi? 



(37) 



Thus, defining C s (w,a s ,l) = C g (u>,a s ), the resummed gluon coefficient function in the massive gluon regularization 
scheme is given by 



C g (u},a s ) 



C A 2 



^ui 2 + 16Cao, s 



1 



Furthermore, from Eqs. (|25[) and (|36f) . we find that 

f = 



1 ^f(r-i) 
o r 



Using Eq. (p3|) . we find that f obeys the DGLAP equation 

d 



r = p DL f, 



din iXp 

where P DL = is the complete DL contribution to all the splitting functions, with 

n 2Cf 

1 o Y 



(38) 



(39) 



(40) 



(41) 



Note that Eq. (|39|) can be obtained by solving Eq. (j40|) with the boundary condition f (w, a s , 1) = 1, and then using 
the projection operator property A 2 = A, i.e. 



/4 

T | u),a a , — 5- 



exp 



in 



P DL (cj : a s )\n 



<j_ 



1 + ^4 exp 



7(w,a s )ln- 



I'f 



(42) 



To check C(uj,a s ) and 7(w, a s ) against full FO results in the literature, we consider the probability distribution g 
for the inclusive production of a gluon in e + e~ annihilation. In this case, Kj = 2Cf, where the factor 2 accounts for 
the fact that the gluons can be emitted either collinearly to the quark of momentum p or the quark of momentum 
p. Up to NNLO, P DL agrees with the FO results for the DL contribution to the splitting functions calculated in the 
literature, as shown in Ref. Q . The NNLO expansion of Eq. 



C g (uj,a s ) 



2C F 
~~Ca~ 



-4CU^+48Cl(^ 



agrees with Ref. [|J. However, it disagrees with the small ui limit of the full NNLO result e.g. in Ref. @, 



C g (ui,a s ) 



2C F 
~Ca~ 



_4CU^+40Ci(^ 



(43) 



(44) 



As already said, this discrepancy is not surprising, given that the result in Ref. Q was obtained using dimensional 
regularization followed by MS factorization, while here the result is obtained using a cutoff to regulate the collinear 
singularities. In section Mil we will calculate the DL contribution in the MS scheme, and show that its NNLO 
expansion agrees with Eq. ([4"4"]). 



III. DOUBLE LOGARITHMS IN THE MS SCHEME 



In this section, we calculate the quantity G, this time using dimensional regularization. This will allow us to 
implement MS factorization. This requires deriving in terms of the variables x a and z a the probability dwi (q a ) in d =/= 4 
dimensions with a zero gluon mass, i.e. Eq. (|A1 1|) with k = q a . Now d d ~ 1 q a = \q a \ d ~ 2 d\q a \siri d ~ 4 6 a dcos9 a dQd-2, 
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where fid is the <i-dimcnsional solid angle. Because of the azimuthal symmetry of Eq. (|A11[) with respect to p, the 
integrand is independent of fld-2i which may therefore be explicitly integrated out using Jdn d = dir d / 2 /T(d/ 2 + 1) 
to give 



dwi(q a ) = 2a s ^' ^ ) \ K I - 2{ - P '. p) \$ a \ l - %t sin' 26 a d\q a \dcos6 a , (45) 
r(l-e) {p-q a ){p-q a ) 

where \x is an arbitrary parameter with the dimension of mass (called the dimensional rcgularization mass) which is 
needed to ensure that the action is dimensionless (note that gfi e is independent of (i and g is dimcnsionlcss). Now, 
performing the change of variables given in Eqs. (|5J) - (|10[) with m g = gives the dimensionally regularized alternative 
to Eq. dTJ), i-e. 

a I \ o {V 2 Y ( 47r ) £ v dx <* dz <* (AR\ 

d Wl (x a ,z a )=2a s ^) _^_ y i f/ _ F3 -_. (46) 

We have omitted a factor (1 — z)~ e = 1 + O(z), since doing so does not affect the DL contribution. In this case, the 
collinear cutoff Z can be set to zero, because the integrals for e < are well-defined even in the limit Z — > 0. Hence, 
Eq. (fT8)) is replaced with 

1 dx\ f 1 dx2 f 1 dxN-i 



xG(x, e) = 5(1 - x)+x-* £ X N I 4^7 / 

JV=1 •' x X l "i 



1 dzi f 1 dz2 f 1 dzN 



l+2c ' - ' / ^.l+2e 

(47) 



x 2 Jxn-2 X N~1 



where 

T(l-c) \Q 

For comparison with Ref. Q, we note from Eq. (|47]l that G(x, e) = G(x, Z = 0, e), where G(x, Z, e) is another quantity 
defined implicitly by the integral master equation 



X = 2C A a s J^^ (£ ) . 



x^G(x, Z, e) = 5(1 -ij+l/ 1 ^ f j%- x'G(x', z', e). 

Jx % J z (z') 1+£ 

Equation (|49p represents the d = 4 — 2e generalization of the four-dimensional master equation given by Eq. (fH)]) 
together with Eq. JT6]). 

In the MS renormalization scheme, the bare coupling a s is related to the renormalizcd coupling a s (fi R ,e) via 
a s S e (n 2 /n 2 R y = a s (/4,e) + 0(a 2 s ), where S t = e < lni ^ E \ so that 

X = 2C A a s (^ e) (j^j [l + (e 2 )] + 0(a 2 s ), (50) 

where the 0(a 2 ) and 0(e 2 ) corrections are neglected because they would only result in terms in our final factorized 
quantities which are sublcading relative to the DLs that we are interested in here. The quantity \x R is a renormalization 
scale, of which X is formally independent. Using the symmetry of the z integrals, we obtain 

f 1 d Zl f 1 dz 2 f 1 dz N 1 ( f 1 dz \ N _ (~1) N 



2 1 Jzi Z 2 Jzx-i Z N 



Similarly, the x integrals are given by 



1 dxi f 1 dx 2 f 1 dx N -i 1 f f 1 dxi \ N 1 (ar 2e - l)*" 1 



Xi J-2 J %N-2 X N-1 



9N{X) ~J X xl+^J XlX l^---J XN _ 2 x^t (N-iy. U x\^) (N-l)\(2e)»-i- 

Therefore, 



1 — (-,, N T- 2e (-r- 2e - "I) 1 *- 1 

G(x, Xl e) =5(l-x) + - £ XN N ! {N ^iy. 2 N- le 2L ■ (53) 
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To calculate the Mcllin transform, wc first define and then evaluate using integration by parts the quantity 



dxx^ix'^ - l) r 



2er 

UJ 



2er 



(J r + I r _i). 



Thus, I r = I r -i/[(jj/(2er) — 1]. Using Iq = 1/ui, we find that 

t-i i 



w rf - 1 



(54) 



(55) 



Using this result, we may evaluate the Mellin transform of the x-dependent part in Eq. (|53[) to obtain 



dxx^ar^ar^ - 1) 



\k-l 



where in the first step we have used Eq. (|54"|) . Thus 
We note that 



(56) 



fe=0 



<;,*,e) = ofi(; i-£; ^ 



(57) 



(58) 



where p F q (a\, ci2, . . . , a p ; b\, 62, • ■ • , b q ; z) is the generalized hypergeometric function. (In fact, qF\ is related to 
the Bessel function I u via 0-F1 l;z 2 /4) = \T(y + 1) / '{z/2) v ]I v {z).) Now, direct factorization of the collinear 

singularities as e — > in Eq. (|57|) by expanding it in e does not seem possible. However, from Eq. (|57j). it is easy to 
check that G satisfies the following simple differential equation: 



G-fG X 



G 



2e 2 '' 



where wc have defined 



(59) 



(60) 



Therefore, knowing from the factorization theorem that the factorization of the collinear singularities in Eq. (|57[) takes 
the form 



G(uj, X,e) = C ^w, a a ([i%,e), ~T'^j ^ ( w ' a sG"|> e), e) , 



(61) 



where the transition function T contains all the e — > collinear singularities, we can solve Eq. (f5T))) to obtain the 
timelike coefficient function C using the fact that it is finite as e — > 0. In the MS factorization scheme, the timelike 
splitting function, 



j(w,a s (fi 2 F ,e)) 



lnr(w,a s (/i|,e),e), 



<iln 

is explicitly independent of e. Using the boundary condition T(uj, 0, e) = 1 for the transition function, we have 

r/4 W«2 



F(uj, a s (fi F ,e), e) = exp 



/ — l{u,a s (q ,e)) 
Jo 1 



Using the equation for the running coupling, 

da s (fi 2 R ,e) 



-ea s {iJ, 2 R , e) - ^n< +2 0«R, e), 



(62) 



(63) 



(64) 



71=0 
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and the fact that a s (0, e) = 0, we can perform a change of integration variable in Eq. (|63[) from q 2 to a s (q 2 , e), which 
gives 



T(uj, a s , e) = cxp 



da. 



To our accuracy, 



T(uj, a s , e) = cxp 



1 f as da. 



(65) 



(66) 



because the inclusion of the terms of 0(a 2 ) and higher in Eq. (|64[) does not affect the DL contribution. For simplicity, 
we will set [Ir = fiF = Q from now on, so that 



and 



X = 2C A a s (Q 2 ,e) 
G(u>, X,e)=C (w, a s (Q 2 , e), e) T ( W) a s (Q 2 , e), e) 



(67) 



(68) 



where C (w, a s , e) = C (uj, a s , 1, e). 

The simplest way to obtain the timelike splitting function 7 and the timelike coefficient function C is to compute 
the left hand side of the differential equation in Eq. (|59p using Eqs. (|6"8]) and (|6"6"]) , and then to compare the two 
inhomogeneous terms. In this way, we will find two solvable implicit equations which will allow us to find closed 
expressions for the DL contributions to 7 and C. So, differentiating Eq. (|68| with respect to X and using the result 



f = --r, 

e 



that follows from Eq. (|66l) , we obtain in the notation of Eq. (|6Cfl) the results 



G 



and 



so that 



G 



C-2-C 



/ 1 UJ /~i 

G 



27- 



c 
c 



+ 7 + 



c 

G' 



(69) 



(70) 



(71) 



(72) 



Note that keeping the factor K e = 1 + 0(e 2 ) in Eq. ([50)1 would result in an O(e ) change in this expression coming 
from the first term of 0(e~ 2 ), thus X e can be safely excluded. Now, comparing the coefficients of e~ 2 and e _1 on the 
right hand side of Eq. (|T2")) with those of Eq. (f5"9"f and noting that 7 is explicitly independent of e gives respectively 



9 Cd 7 X 



and 



dlnC 

97 



1 1 



2 7 + 



(73) 



(74) 



Equation ()73j) is the same quadratic equation for 7 as the one appearing in Ref. [H, [l4j], whose unique solution 
compatible with perturbation theory is Eq. (f32j) . Solving Eq. (|74| for C gives 



C 



My) 



(75) 
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where A{uj) is an unknown constant of integration. Using Eq. (|32l) for 7 and the condition C(uj, 0) 
perturbation theory, we obtain the result A{oj) = t/uj/2, so that, finally, we find 



From Eq. (|37p . we obtain 



C(uj,a s ) = 



Ki 



1 



1 + 16CV 



1 from 



(76) 



C g (u>, a s ) = — 



This is the main result of this paper. As us — > 0, C g is finite and is approximately given by 



C g (uj,a s ) 



Ki 



2C A {C A a s y 



To obtain C„ in x space, we expand it in a s to obtain 



C g (u>,a s ) 



El 

Ca 



16C A a, 



- 1 



Kj_^ (_l)r (i/ 4)r fiQC A a £ 
C A ^ 



r=l 



\ 



where (a)k = T(a + k)/T(a) are the Pochhammer symbols. Then, using the result in Eq. (|24[) . we obtain 



C g (x,a s ) 



K T 1 



(l/4) r 



16CUa s In — 

a; 



C A! rlni^(2r-l)! r! 
Identifying the sum gives the result for the DL contribution to the gluon coefficient function in x space, 



C g (x, a s ) = -AKids 



In - 



1^2 



5 3 

^; 2' 2; -4C*Aa s ln 2 x 



(77) 



(78) 



(79) 



(80) 



(81) 



Finally, a gluon coefficient function which is valid from both large to small x and formally consistent with the FO 
and SGL resummed approaches is the DL-resummed NLO gluon coefficient function, given by 



C, DL + NLO (x, a.) = C^° (x, a s ) + C^ (x, a.) - AKia ° lnx , 



(82) 



where C® 1 " is given by Eq. (|5Tj) . Cg^ ho by the NLO expansion of the FO result for C g (e.g. for e + e _ annihilation, see 
Eq. (|84|) below), and the last term in Eq. (|82|) prevents double counting of the DL at NLO. 
In the case of e + e~ annihilation, we have Kj = 2Cf, so that 



c. 



1 



{i + i&c A ^Y 



(83) 



In contrast to Eq. (|38p . Eq. ([53^1 now agrees with the FO results in the literature, i.e. it is equal to Eq. (|44[) when 
expanded to NNLO. More generally, Eq. (|77| resums all the DLs of the gluon timelikc coefficient function for e + e~ 
annihilation in the MS scheme. The gluon timelike coefficient function only starts to contribute at NLO in QCD. The 
contribution at NLO in x space is given by [TH, [l|| 



(84) 



C„ NLO (x, a s ) = a s C F \ 2^-t-^ — ^- [ln(l - x) + 2 In a?] 



In Fig. [TJ we plot the NLO gluon coefficient function in Eq. (|84|) . together with the pure resummed result in Eq. (|8Tj) 
and the combined DL+NLO result in Eq. (|52"j) as functions of ^ = ln(l/a;). Note that the gluon coefficient function 
is negative. We see that the DL+NLO result of Eq. (|82|) interpolates well between its 0(a s ) FO result (Eq. (|84|l) at 
large x and the pure DL contribution (Eq. (|80l) ) at small x. 

We note that Eq. (fB"2")) does not show a DL divergent behaviour, or in fact any kind of divergent behaviour, for 
small x, in contrast to C^ LO , and this fact is illustrated in Fig. [T] in the large £ region. In fact, the quark coefficient 
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1 2 3 4 5 

FIG. 1: The gluon coefficient function in e + e — annihilation 
at NLO in Eq. 1841 1 (blue lower line), the complete DL con- 
tribution to the gluon coefficient function in Eq. I|80| l (violet 
line) and the DL-resummed NLO result in Eq. 18211 (orange 
line). The value a s = 0.18 was used, because this corresponds 
to the result for a s (Q 2 ) at LO using Q = 14 GeV, the value 
for Q used in Ref. Q]. 



-xg h {x,Q 2 ) 




1 2 3 4 5 



FIG. 2: Contribution to the inclusive single hadron produc- 
tion cross section in e + e _ annihilation due to the gluon chan- 
nel corrections in Eq. I|85j l at NLO and with the inclusion 
of the DLs at all orders in perturbation theory. As in Fig. 
Q] the value a s (Q 2 ) = 0.18 was chosen in the calculation of 
C g {z, a s (Q 2 )), and the gluon FF was taken to be that obtained 
at Q = 14 GeV in the LO global fit of Ref. Q]. 



function at NLO is also free of DLs and SLs (see e.g. Refs. [HI, (HI). In other words, there are no large logarithms at 
small x at the NLO level after rcsummation of the DLs in the gluon coefficient function. Thus it is possible that the 
LO analysis done in Refs. [2, El ma y now be extended to NLO, and for a similar range of data. In Fig. ([2]), we have 
estimated the contribution to the inclusive single hadron production cross section coming from the gluon coefficient 
function in the MS scheme, which is calculated according to the a;-space convolution formula, 

«"<^> = ik'^^ " II ^(s<, (ftf) ■ (85) 

where all scales have been chosen equal and Dg is the gluon FF. The quantity (1 / 'da B)(da h / 'dx) has been evaluated 
with C g equal to the NLO result in Eq. (f54")l , and to the improved DL+NLO result in Eq. (|5^|) . which includes DLs 
at all orders. The gluon FF has been chosen to be that obtained at yfs = 14 GeV in the global LO fit in Ref. Q, in 
which all the DLs in the anomalous dimension were resummcd. As expected, this contribution is negative. Even if 
the correction seems to be large, this does not imply a correction of the same amount in a NLO global fit, because 
the gluon FF obtained in a LO analysis carries a large theoretical error. Nevertheless, at the same time, this suggests 
that a sensible negative correction around and above the peak for the total hadron multiplicity is to be expected, thus 
necessitating a complete analysis using also a DL-resummed gluon coefficient function [171 ] . 

Finallv, we note that, given that the SLs in the anomalous dimensions calculated in the MG scheme are already 
known [18[, all that is needed to calculate the SLs in the anomalous dimensions calculated in the MS scheme is the 
relation between these two schemes at the DL level. This relation is partially constrained by the fact that C g is now 
known in both schemes (Eqs. (|38|) and ([77)) ). Thus, although a complete determination of the SLs in the anomalous 
dimensions calculated in the MS scheme is still not possible, some constraints may now be imposed on them which 
e.g. could be tested against FO results. Further studies will be left to a future publication 



IV. CONCLUSIONS 



In this paper, we have computed and resummed all the DLs present in the gluon coefficient function. To resolve 
a mismatch between existing resummation formulae in the literature [3| and recent higher order corrections up to 
0(a 2 s ) (5T4l0j. we have obtained the DL resummed result by applying dimensional regularization and then subtracting 
the collinear singularities according to the MS factorization scheme. The FO expansion of our new result agrees with 
the DLs of the full NNLO calculation. We have seen that the resummation of the DLs in the anomalous dimensions 
and coefficient functions is sufficient to eliminate all large logarithms at NLO that could destroy perturbation theory 
at very small x values. This could enable us to obtain observable quantities at NLO from data down to the same 
small values of x as were previously treated at LO in Refs. 0, Q. We expect the impact of the DL-resummed gluon 
coefficient function to be important, but not necessarily to modify drastically other results previously obtained, and 
we have pointed out the necessity of a complete DL+NLO analysis [l7| . 
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Appendix A: Factorization of multiple gluon emission 

Here we briefly sketch how the gluon emission probability factorizes in the kinematic region from which both soft 
and collinear DLs arise. ("Factorize" in this appendix refers to the vanishing of colour correlations.) Here we recall 
only the features which are needed to understand the derivation of Eqs. ([TT]) and (|4*o]) for the gluon emission probability 
in d = 4 and d = 4 — 2e dimensions, respectively, with e < 0. For a more general and complete treatment about this 
factorization for infrared sensitive quantities in QCD jets, we refer the reader to Ref. [l9|, of which the following is a 
summary. 

We work in d = 4 — 2e dimensions. We start by recalling that the QCD eikonal amplitude £i(p,e(k)) associated 
with a soft gluon with momentum k and polarization vector e M (fc) that is emitted from a parton I (I = Q, if the 
parton is a quark, Q if it is an antiquark, and g if it is a gluon) of momentum p is given in the case where |fc| <C \p\ 
by the result 

£ I {p,k)=g^T I j^(p,k)e^(k), (Al) 

where jfi(p,k) = Pn/(jp ■ k). Tj are the colour charge matrices which act on colour space only, and their matrix 
elements for the emission of a gluon of colour index a are given by (Tj)^- = tfj (the colour matrix in the fundamental 
representation) if I = Q, by (Tj-)|y = —VL if I = Q and by (T/)£ c = ~if a bc (the colour matrix in the adjoint 
representation) if I = g. We remind that the colour-charge algebra is given by: 

Tj • Tj = T, 7 ■ Tr if I ^ J; T? = Kj 1, (A2) 

where K q = Kg = Cf = {N% — l)/2N c , K g — Ca = N c and 1 is the identity colour charge operator in the 
fundamental and adjoint representations if / = Q, Q and / = g, respectively. If we have many gluons in the final 
state, then radiation of gluons by other gluons can occur, and the approximation in Eq. (|A1[) can be applied iteratively 
when the momenta of these gluons are strongly ordered. 

We now want to consider the more general physical case in which there are N soft gluons emitted with momenta 
k a , where a = 1, 2, . . . , N. We make the assumption that \k\ \ <?C \k2\ -C ■ • • <C |fcjv| <C |p|. Because we restrict our 
phase space to only one strongly ordered region, we treat each of the N gluons as distinct when summing over the 
modulus squared amplitudes. Because the \k a \ are negligible with respect to \p\, we also have to add a recoiling parton 
having momentum p with p w — p in order to respect momentum conservation. Because the gluon with momentum 
k\ is the softest gluon, the emission amplitude £t ot (p,p, {k\, &2, . . . , k;sr}) for the whole process can be calculated by 
calculating the emission amplitude £tot{p,P, {fe, ^3, ■ • ■ , fc/v}) and then attaching the gluon with momentum k\ and 
polarization e(ki) to each of the final state real partons. The result is 

N+2 

£tot(p,P,{ki,k2,---,k N }) = ^2 £i(a)(k a ,ki), (A3) 

a=2 

where 1(a) = g for a = 1, 2, . . .N, I(N + 1) = Q, I(N + 2) = Q and where we define k^+i = p and kx+2 = P- Now, 
performing the square modulus of this amplitude and summing over the polarizations of the final state gluons, we 
determine the insertion operator for the soft gluon to be 

N+2 

I(p,p,{fc 1 ,fc 2 ,...,fc JV }) = ^|f tot (p,p,{fc Q })| 2 = 3 V 2e E T;wT ;(fl j"(* 0) ii)j , '^ l fc 1 )dW(ti), (A4) 

£ a, 13=2 

where dju) is the on-shell polarization sum which, when parametrized by a gauge null vector n M , reads 

4" ) (fc) = -.^+ fc ^ + , fc '^ . (A5) 
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Substituting this into Eq. (|A4[) . we find that the insertion operator can be written as 

I(p,p, {h,k 2 , • ■ • , k N }) = -g 2 [i 2e J M (p,p, {h, fc 2 , . . . , k N })3^(p,p, {ki, fc 2 , . . . , k N }), (A6) 



where 



N+2 



/Iff 1 n 11 \ 

J»(p,p,{k 1 ,k 2 ,...,k N }) = Y,Ti{«) [k^-^TY)- ( A7 ) 

a — 2 ^ ' 

We note that the quantity J M (p,p, {fci, k 2 , ■ ■ ■ , fcjv}) is gauge invariant if we assume for the final state of the process 
without the gluon with momentum fci (described by £tot(p,.P, {k 2l k 3 , ■ • ■ , k^})) the color singlet condition 

N+2 

E T '(«) = °> ( Ag ) 

a=2 

which implies that the gauge vector n M can be chosen in a way that is convenient for calculations. Due to colour 
correlations among the emitted gluons in Eq. (|A4|) . in general we expect that the emission probability does not 
factorize. Nevertheless, this does happen in the angular configuration from which the DL contribution arises, namely 
the configuration in which fc/v is approximately either parallel to p or anti-parallel to p (i.e. parallel to p), and the 
gluons of momentum k a for a = l,2,...,iV — 1 have much smaller angles of emission relative to the angle of emission 
of the gluon of momentum fcjv, and are approximately either all parallel or all anti-parallel to fc^r. Because they then 
form a jet around cither p and p, i.e. either fc^ ~ KaP 11 or fc^ ~ HaP 11 for a = 1,2, ...,N— 1, where n a are some 
parameters that are independent of the spacetime index /i. For definiteness, we will assume that fcjv is approximately 
parallel to p and then later multiply by a factor 2 to account for the case that kjsi is approximately parallel to p. 
In this configuration, we find that the gluon with momentum k\ factorizes and carries the total colour charge of the 
jet: We choose n M = p in Eq. (|A7[) for convenience, note that k^/k a ■ k\ = p^/p ■ fci for a = 2, 3, . . . , N + 1, and 
then use the gauge invariance condition in Eq. (|A8[) to write y^ Ja —2 = — T/(jv+2)- If the initial state partons of 

momentum p and p are gluons, then Ti(n+2) — T g . If they are a quark and an antiquark, then, assuming that there 
are gluons already attached before the gluon of momentum k\ is attached, these quarks must form a colour octet so 
that Ti[n+2) = T g again. Therefore, using the second result in Eq. (IA2[) . we find that the insertion operator in Eq. 
[6|) now simply becomes 



2{P-P) 
(p- ki)(p- fei)' 



I(p,p, {k u k 2 , . . . , k N }) = 5 2 M 2e K g j- u - . (A9) 



Hence the single gluon emission probability for this angular configuration factorizes according to 

da N (p, p, fci , k 2 , . . . , k N ) = da N ^x(p, p,k 2 ,k 3 ,..., k N )dw g (fci) , (A10) 



where 

dwj(k) = g 1 ^ Kj , j- n (ah) 



,2,>. 2(p-p) d d -'k 



9 [l 1 {p-k)(p-k) (2 7 r)' i - 1 2fc 

is the probability of emission of a single soft gluon of momentum k from a parton /. Here, da^ is the A^-gluon 
emission cross section and the last factor is the usual d-dimensional phase space measure. In Eq. (|A10|) . all momenta 
k a arc fixed, where a = 1, 2, . . . , N. Note that the gluon with momentum k± is distinguished from the others only by 
the fact that it is the softest. Thus, because the gluon with momentum k 2 is the softest one in the process described 
by daN-i(p,p, k 2 , . . . , fcjv), by the same reasoning we have 

rfo-Ar_i(p,p, k 2 ,k 3 ,...,k N ) = d<7 N - 2 (p,p,k 3 ,k4, . . . ,k N )dw g (k 2 ). (A12) 

Therefore, iterating this procedure, we find 

da N (p,p, fci, k 2 , . . . , k N ) = dcr B {p,p)dwg(k 1 )dwg(k 2 ) . . . dw g (fcjv-i)(iw/(Ar + 2)(fciv)- (A13) 

Note that, in the last step, we have used dai(p,p,kpf) = doB(p 1 p)dwi(N+2){kN), where dwi^M+2) dw g necessarily 
because Tj^ N+2 ^ = T q or Tq if the initial state partons are a quark and an antiquark: they must be in a colour singlet 
state due to the fact that there are no gluons in the final state before the gluon with momentum k^ is attached. 
Note that the same result in Eq. (|A13I) would have been obtained if we had used a different strong ordering for the 



14 



momenta k a for a = 1, 2, . . . , N — 1. This is a special case of the invariance under choice of strong ordering that is 
proved in Ref. (Tlj . 
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